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I.  Introduction 


Recent  work1  J  in  a  variety  of  areas  of  surface  spectroscopy  indicates 
that  it  is  possible  to  modify  in  a  radical  manner  the  spectroscopic  signals  of 
a  molecule  by  placing  it  near  a  solid  surface. 

Such  modifications  can  be  caused  by  three  electromagnetic  effects. 

(1)  The  surface  is  polarized  by  the  incident  beam  and  the  molecule  interacts 
with  both  the  incident  and  the  polarization  field.  The  local  light  intensity  is 
thus  higher  and  the  probability  of  a  molecular  excitation  is  enhanced. 

On  a  flat  surface  this  mechanism  produces,  at  best,  a  "minor  enhance- 

m 

ment"1  J  similar  to  that  obtained  by  placing  the  molecule  near  a  mirror 
(whose  quality  depends  on  the  solid  material  and  the  frequency  of  the  incident 
light).  For  curved  surfaces,  however,  the  light  can  excite  an  electromagnetic 
resonanc  Jtf.g] 

in  the  solid  which  has  a  strong  evanescent  field  near  the  surface. 
This  causes  a  substantial  enhancement  of  the  local  intensity  and  of  the  prob¬ 
ability  of  molecular  excitation.  Furthermore,  the  curvature  of  the  surface 
enhances  the  field  even  in  the  absence  of  a  resonance. 

T2l 

(2)  The  same  mirror  effect  can  produce  a  minor  enhancement1  J  of  the 
emission  from  an  excited  molecule.  The  process  is  more  interesting,  however, 
if  the  molecular  emission  frequency  coincides  with  that  of  an  electromagnetic 
resonance  of  the  surface.  The  oscillating  dipole,  representing  the  emitting 
molecule,  can  drive  and  excite  the  resonance.  The  spectroscopic  outcome  of 
this  depends  on  whether  the  resonance  prefers  to  radiate  or  to  use  its  energy 

to  heat  up  the  surface.  If  the  former  is  the  case  the  total  emission  is  sub¬ 
stantially  brighter. 

(3)  The  spectroscopic  properties  of  the  molecule  are  also  affected  by 


the  molecule -surface  energy  transfer.  This  diminishes  the  lifetime  of  the 
molecule  and  it  is  thus  an  "undesirable"  effect,  unless  fluorescence  quenching 
is  of  interest.  Mathematically  the  rate  of  transfer  is  related  to  the  "image 


field"  of  the  molecular  dipole  near  the  surface,  ^2’  10^  which  takes  into 
account  the  fact  that  the  dipole  polarizes  the  surface  and  oscillates  under  the 
influence  of  the  polarization  field. hi  doing  so  it  must  use  a  part  of  its 
power  to  sustain  the  polarization  of  the  solid  and  this  is  taken  from  the  power 
that  would  have  been  radiated  if  the  molecule  was  in  vacuum. 

It  is  likely^- that  other  effects  (i.e.  "molecular"  effects,  related  to 
chemisorption)  are  also  at  work  in  surface  enhanced  spectroscopy.  However, 
we  concentrate  here  on  the  role  of  the  electromagnetic  effects.  These  are 
controlled  by  a  number  of  important  factors. 

(a)  The  shape  of  the  surface  determines  the  frequency  of  the  resonance 
and  the  magnitude  of  the  evanescent  (i.e.  near)  field.  This  can  be  excited  by 
using  an  ATR  (Attenuated  Total  Reflection)  configuration,^  gratings,  ^ 
spheres  and  ellipsoids,  ^  metal  islands,  two  spheres, small  random 
roughness ^  and  a  "coarse"  rough  surface. 

(b)  High  evanescent  fields  can  be  sustained  only  by  materials  which 
are  unable  to  damp  the  electromagnetic  resonance.  The  ones  found  useful 
so  far  are  Au  and  Cu  (~red),  Ag  (blue  and  UV),  SiC  and  InSb  (infrared). 

One  of  the  recent  direction,  taken  in  this  area  is  a  search  for  new 
surface  geometries,  which  can  permit  higher  enhancement  factors,  greater 
frequency  tunability  and  can  function  with  a  broader  range  of  materials. 

From  the  point  of  view  of  the  theorist  one  would  prefer  systems  which  can  be 
prepared  in  the  laboratory  and  also  permit'an  accurate  computation  of  the 
electromagnetic  effects,  in  such  cases  the  difference  between  the  experimental 
observations  and  the  predictions  of  the  electromagnetic  model,  if  any,  can  be 
attributed  to  "molecular"  effects. 


In  what  follows  we  propose  such  a  geometry.  This  consists  of  a 
perfectly  conducting  sphere  "suspended"  above  a  flat  surface  which  has 
an  arbitrary  dielectric  constant.  This  system  can  be  prepared^  by 
directing  on  the  flat  surface,  in  a  vacuum  system,  a  molecular  beam  of 
microscopic  metal  spheres  (radius  ~300A  or  less)  in  a  carrier  gas,  and  a 
beam  containing  the  molecules  on  which  one  would  like  to  do  enhanced 
spectroscopy.  The  temperature  should  be  adjusted  so  that  the  molecules 
of  interest  will  freeze  on  the  surface,  forming  a  matrix  which  traps  the 
metallic  spheres.  Thus  we  obtain  a  flat  surface,  covered  by  a  molecular 
matrix  containing  perfectly  conducting  spheres.  This  structure  can  be 
prepared  for  a  variety  of  purposes.  One  might  use  a  carrier  gas  which 
adsorbs  on  the  surface  of  the  spheres  and  trap  the  spheres  in  a  spectro¬ 
scopically  inert  matrix,  on  a  flat  surface  that  can  provide  the  enhancement; 
in  such  a  configuration  we  can  study  the  spectroscopic  properties  of  the 
molecules  adsorbed  on  the  spheres.  Another  possibility  is  to  preadsorb 
molecules  on  the  flat  surface  and  then  use  spectroscopically  inert  carrier 
gas  to  form  a  matrix  which  entraps  the  spheres  near  the  flat  surface.  This 
can  be  used  to  study  the  spectra  of  the  molecules  adsorbed  on  the  flat  surface. 

The  mechanism  for  enhancement  in  this  geometry  can  be  qualitatively 
understood  by  the  description  of  the  following  sequence  of  events.  Assume 
first  that  the  sphere  is  absent.  The  laser  field  cannot  excite  the  surface 
plasmon  of  the  flat  interface  because,  if  it  were  to  do  so  it  would  violate 
parallel  momentum  conservation.^®^  The  presence  of  the  sphere  breaks 
the  translational  symmetry  along  the  surface  and  parallel  momentum  con¬ 
servation  ceases  to  be  a  kinematic  requirement.  The  plasmon  can  now  be 
excited.  The  mechanism  of  this  excitation  can  be  understood  in  the  following  manner. 


The  incident  light  polarizes  the  sphere  and,  for  large  sphere  surface 
distances,  the  polarization  field  acting  on  the  surface  is  that  of  a  dipole. 

This  is  known^^  to  excite  the  plasmon.  If  the  sphere  is  brought  closer 
to  the  surface,  it  is  polarized  by  both  the  incident  laser  field  and  by  its 
own  image  field.  The  latter  is  spatially  inhomogeneous  and  can  excite  the 
quadrupole  in  the  sphere,  which  in  turn  can  also  excite  the  plasmon.  The 
plasmon  interacts  with  the  dipole  and  the  quadrupole  and  its  frequency  is 
shifted.  Thus  the  properties  of  the  resonances  of  this  structure  are  dif¬ 
ferent  from  those  of  the  plasmon  in  the  absence  of  the  sphere.  If  the 
sphere  moves  closer  to  the  surface  the  octupole,  etc.  also  play  a  role. 

The  outcome  is  that  the  excited  plasmon  interacts  with  the  polarized  sphere 
and  builds  "new''  surface  resonances  which  enhance  the  electric  field 
between  the  sphere  and  the  surface  and  whose  excitation  frequencies  are 
very  different  from  those  of  the  surface  plasmon.  This  large  local  field 
can  be  used  to  carry  out  surface  enhanced  spectroscopy. 

The  calculation  of  the  local  field  is  rather  straightforward  in  the  limit 
when  the  distances  involved  (the  radius  of  the  sphere  and  the  plane-sphere 
separation)  are  much  smaller  than  the  wavelength  of  light.  When  this  is 
the  case  retardation  can  be  neglected  and  the  scalar  potential  can  be 
calculated  by  solving  Laplace's  equation.  The  solution  is  obtained  by  using 
bispherical  coordinates  and  the  method  of  separation  of  variables,  as  exempli¬ 
fied  for  many  cases  in  the  book  of  Morse  and  FeshbachJ*^ 

The  remainder  of  the  paper  is  organized  as  follows.  In  the  next  section 
we  specify  the  geometry  and  outline  the  procedure  used  to  solve  the  Laplace 
equation.  Section  III  presents  numerical  results  chosen  to  illustrate  the 
properties  of  the  electromagnetic  resonances  of  this  particular  structure.  In 
the  last  section  we  comment  on  possible  applications  and  compare  this  system 
to  other  candidates  for  surface  enhanced  spectroscopy. 


II  The  Model 


n.a.  The  geometry  and  the  materials 

The  geometry  considered  here  is  shown  in  Fig.  1.  It  consists  of  two 
half  ^spaces  separated  by  a  planar  surface.  The  "optically  active"  half-space 
1,  of  dielectric  constant  c^(oj),  earn  sustain  a  surface  excitation  (i.e.  Re  €(co) 
can  become  negative  in  some  frequency  range),  such  as  a  surface  plasmon, 
polariton  or  exciton.  The  other  half-space, in  which  the  sphere  is  imbedded, 
is  "optically  inert":  its  dielectric  constant  is  real  and  frequency  inde¬ 
pendent.  The  sphere  is  located  inside  the  inert  medium  at  a  distance  D 
from  the  flat  surface.  It  has  a  large  dielectric  constant  (i.e.  good  conductor) 
and  its  radius  is  RQ. 

This  structure  is  driven  by  an  external  electromagnetic  field  of  frequency 
C t} .  The  wave  vector  5  is  contained  in  the  x-z  plane  and  makes  an  angle  0Q 
with  the  normal  to  the  surface  (i.e.  the  z-axis).  To  describe  the  polarization 
we  use  the  nomenclature  associated  with  flat  interfaces:  the  p-polarized 
beam  has  the  electric  vector  in  the  plane  of  incidence;  the  s-polarized  beam 
has  the  electric  vector  perpendicular  to  the  plane  of  incidence  (parallel  to 
the  surface). 

n.b.  The  method  of  solution  of  the  Laplace  Equation. 

We  compute  the  fields  in  the  inert  medium  for  the  case  when  RQ  and 

D  are  smaller  than  the  wavelength  2irc/u}  (Rayleigh  limit).  In  this  case  we 

2 

need  only  to  solve  the  Laplace  equation  7  <5=  0  for  the  scalar  potential. 


The  expression  for  the  potential  is  chosen  to  be  of  the  form: 


$  =  $ext  +$*  +  <J>3  .  (1) 

Here  $ext  corresponds  to  the  laser  field  and  $  to  the  (Fresnel)  field 
reflected  by  the  flat  surface.  4*s  is  unknown  and  takes  into  account  the 
modifications  introduced  by  the  presence  of  the  sphere. 

The  external  field  (specified  in  Fig.  1)  is  given  by 


.ext,-* . 

$  (r  )  =  - 


£?i  -*  -iU)t 

E  •  r  e 


(2) 


The  amplitude  E1  of  the  electric  field  is  taken  to  be  independent  of  the 
spatial  coordinate  since  this  is  roughly  true  over  the  region  of  interest 
(around  the  sphere).  We  have 


and 


E1  =  E1  =  0;  E1  =  E„ 
x  z  y  o 

(for  s -polarization) 


E*  =  -  E  cos  9 .  E*  =  E  sin  0  ,  E*  =  0, 
x  o  c t  z  o  o  y 

(for  p-polarization). 


(3) 


(4) 


where  EQ  is  independent  of  the  spatial  coordinates . 

The  potential  $ ,  caused  by  the  presence  of  the  flat  surface,  is  chosen 

r  i2l 

to  yield  the  electric  field  given  by  Fresnel  formulae, 1  J 


f  -.r  _ 

$  =  -  E  •  r  ,  for  z  >  0  , 

and 

<I>f  =  -  E1  •  r,  for  z  <  0  , 


(5) 

(6) 


r 

Since  only  the  reflected  amplitude  E?r  appears  explicitly  in  the  present 

calculation  we  reproduce  the  formulae  here.  For  s -polarized  light: 

Er  =  Er  =  0 
x  z 

and 

Ey  s  Eo  (cos  9o"  ^cos  0o+  s  Eors 

with 

a=  "  sin2  0Q 

For  p-polarized  light: 

=  Eq  cos  80  {(c^)  cos  0O  -  a)\  {(ej/e^cos  0Q+a^  s 

(8) 

hEqCOs  0q  rp  , 

Ez  =  K  tan  9o  * 
and 

Er  =  0. 

y 

-•  r  i2i 

The  transmitted  field  E*  is  also  given  by  Fresnel  formulae,  J  and 

we  need  it  when  we  impose  the  boundary  conditions. 

The  presence  of  the  sphere  affects  the  total  field  and  these  modifications 

s  s 

are  described  by  $  .  To  write  an  adequate  form  for  the  unknown  $  we  use 

the  bispherical  coordinates^ ^  n,  rj,  cp  .  The  x-y  plane  (i.e.  the  flat  surface) 

was  chosen  to  coincide  with  fj,  =  0,  while  the  surface  of  the  sphere  coincides 

with  n  -  nQ  (nQ  >0).  The  values  of  nQ  and  the  scale  factor  cQ  are 

determined  by  RQ  and  D  through 

M0  =  arccosh  (l+D/RQ)=fn  |[l+(D/Ro)]  +  [2(D/Rq)+(D/Ro)2]1/2}  (9) 

c 


{(1  +  D/Rc)2  -  l} 


1/2 


(10) 


L 


g 

The  expression  for  $  in  the  inert  medium  (0  <fx<^Q)  has  the  form 

oo  +co 

^  jA^exp  [(n+l/2)|i]+ 
n$^|  m=-oo  (U) 

exp  [-(n+l/2)ju]|  Y™(cos  r\,  (ft)  . 

Inside  the  active  medium  (/ j.  <  0)  we  write: 

oo  +oo 

$S(M,+7,<p)  =  F  D^exp  [(n+1/2)^]  Y™(cos  77,p)  (12) 

n  s  |m  J  m=-oo 

The  definition  used  for  the  spherical  harmonics  Y^1  is  that  of  Morse  and 

FeshbaclJ-*^  and  F  =  (cosh  u-  cos 

We  have  chosen  these  functional  forms  for  the  potentials  because  they 

are  general  solutions  of  the  Laplace  equation  with  the  correct  behaviour  at 

infinity.  The  unknown  coefficients  Am ,  Bm  and  Dm  are  determined  by 

n  n  n 

imposing  the  boundary  conditions 

$(M=  0+,tj,  <p)  =  $(M=  0",  rj,  (p),  (13) 

CQ  |f-  (M=  0+,  77,  (ft)  =  ?1(co)  -||-  (4=<T,  77.  (ft)  (H) 

and 

$  (WQ,  7),  (ft)  =  VQ  .  (15) 

The  first  two  match  the  potential  and  its  derivative  at  the  flat  interface. 

The  inert  half  space  is  at  n  =  0+  and  the  active  one  at  fi  =  0~.  The  third 
condition  recognizes  the  fact  that  the  field  inside  the  perfectly  conducting 
sphere  is  zero  and  that  the  surface  of  the  sphere  is  a  surface  of  equal 
potential  VQ.  The  potential  VQ  is  determined  from  the  requirement  that 
the  total  charge  on  the  surface  of  the  sphere  is  zero  . 


II. c.  The  results  for  $ 


Implementing  this  procedure  yields  the  following  results  for  the 

coefficients  Am,  Bm  and  Dm . 

n  n  n 

For  s -polarized  light 

A?=Bn=Dn=°  am^l 

(16) 

An  =  A^  =  i(1+rs>coEo/2  |exP  [ (2n+l ) MQ]  +  X ^ 

( 

-1 

•  |47rn(n+l)/(2n+l)|  . 

(17) 

with 

T>  1  =  "  fi<“>)]  [co  +  Cjta)]"1. 

(18) 

Furthermore 

B* 1  =  X  i  A*  ^ 

n  In 

(19) 

and 

(20) 

For  p-polarized  light  we  obtain 

A°=  [47T/(2n+l)]1/2  {^Vo+  coEo/T(l+r  )(2n+l)sin  6  ] 

Xr  w  / 

r 

’  {exp  [(2n+l)^]  +3?!!”1  , 

(21) 

An=  -An=/2  coEo[4Jrn(n+1)/(2n+1)]1/2(l-rp)  cos  eQ 

.  {exp  [(2n+l)  Wq]  +X  jJ.  , 

(22) 

Bj*1  -^A®-*1. 

(23) 

=  (i  +jrl)A®-« . 

(24) 

and 


aI£  =  bT=dT=0  U  °*±1 


(25) 


The  quantity  VQ  is 


V  -'oV'V  si“  8„  [e*p((2n+l)(lo)+-t1] 


-II 


(26) 


a=l 


n.d.  The  resonances  of  the  structure. 

An  examination  of  the  equations  giving  the  coefficients  A^1,  B^1 ,  and 
D^1,  for  both  s-  or  p-polarized  light,  shows  that  the  potential  $>s  might 
have  special  properties  (both  inside  and  outside  of  medium  1)  at  frequencies 
0^  satisfying 

(e  (co  )-l\ 

n  =  0,1,2, ...  .  (27) 

When  this  condition  is  fulfilled  the  quantity  { exp [(2n+l)^tQ]  (to)}-1  , 

appearing  in  the  coefficients  A^1,  B^1,  and  D™,  takes  its  largest  value, 

{  1  .  If  Im  (u^)  is  small,  the  potential  $s  can  be  sub¬ 

stantially  enhanced  at  the  frequencies  con» 

It  is  instructive  to  examine  the  frequencies  a)n  for  the  case 
when  the  dielectric  constant  of  the  active  medium  is  given  by  a  Drude  model 
lacking  dissipation,^-14^  for  which  (w)  =  1  -  (u^/co)2.  Here  tOp  is  the 
plasmon  frequency.  For  this  simple  case  Eq.  (27)  gives 

wn  =  —E-  jl  -  exp  [(2n+l)Mo]  .  (28) 

/  2 


If  (2n+l)^0>>l  the  resonance  frequency  “*  t*Jp/ which  is  the 
frequency  of  the  surface  plasmon.  However,  in  the  other  extreme,  when 
(2n+l)/i0  is  very  small,  the  resonance  frequencies  are 


.  . 

n  /T 


{(2n+l)Mo|] 


The  lowest  resonance  frequency,  for  low  (2n+l)^,is  at  <^-0=  • 

From  Eq.  (9)  we  see  that  this  limiting  case  appears  when  X  =  D/Rq<<1 
and  gives  the  lowest  frequency  0)n=0=  Ct^(\/2)  .  If  the  sphere  touches  the 

surface^  con=0  ■*  0. 

For  real  materials  the  analysis  is  more  complicated  and  depends  on 
the  dielectric  constant  Cj (to)  of  the  material  and  the  properties  of  its 
surface  resonance.  A  discussion  of  a  few  practical  cases  can  be  found  in 
Section  HI.  We  sho^y  there  that  even  though  the  coefficients  An,  Bn  and  Dn 
are  maximized  when  o)=wn,  this  does  not  necessarily  lead  to  an  enhanced 
value  of  the  potential. 

n.e.  The  electric  fields  and  the  local  intensity. 

The  electric  field  in  the  space  occupied  by  the  inert  medium  is 
computed  from 


E.  =  -  ~  i  =  u,  77,  <p,  0  <n  <  nQ  ,  (30) 

witht11!  h^  =  h^  =  cQ  [cosh  u  -  cos  77]” 1  and  h^=  cQ  sinr?[cosh  /i-cos  T)]-1. 

K  we  place  an  isotropic  molecule  in  this  field, a  dipole  p(oo)  =  ttE  (ui) 
is  induced  and  this  leads  to  a  spectroscopic  signal  from  the  molecule. 

In  most  cases  the  intensity  measured  in  a  spectroscopic  experiment  is 
proportional  to  p(o>)  •  p(u))*~aa*  £(<*>)•  S  (cc).  The  quantity  E  •  E  is 
thus  a  key  electromagnetic  quantity.  Since  we  have  many  parameters  (such  as 
ou  80.  the’  polarization,  the  point  where  the  field  is  calculated,D/Ro,  Cj  and 
€q)  we  prefer  to  plot  this  "local  field  intensity": 


(31) 


rather  than  give  information  concerning  the  direction  of  E,  which  is  not 
qualitatively  important  for  our  present  purposes. 


H.f.  The  approximations  involved. 

The  use  of  an  electrostatic  calculation  involves  certain  approxima¬ 
tions  which  we  want  to  state  explicitly. 

(1)  The  use  of  the  Laplace  equation  to  compute  the  electric  field 
is  normally  justified  if  the  linear  dimensions  of  the  solid  are  smaller 

than  the  wavelength  of  the  incident  light.  In  our  case  the  flat  interface  between 

the  inert  and  the  active  media  does  not  satisfy  such  a  requirement.  For 

this  reason  we  have  used  the  equation(l)  which  contains  explicitly  (through 

the  exact  solution  of  Maxwell's  equations,  for  the  case  when  the  sphere 

is  absent.  We  use  an  electrostatic  (Rayleigh)  approximation  only  for  the 

computation  of  $  .  This  is  the  potential  caused  by  the  presence  of  the 

the 

sphere.  It  contains  the  effect  of  the  polarization  of  the  sphere  andA"image 
field"  induced  by  the  polarization  of  the  sphere.  The  latter  is  necessary 
because  the  flat  surface  is  present.  This  polarization  is  confined  to  a 
region  whose  linear  extent  is  of  the  order  2(2RQ  +  D).  When  this  is  smaller 
than  2 ffc/co  we  expect  the  electrostatic  approximation  to  hold. 

(2)  The  assumption  that  the  radiation  caused  by  the  polarization  of  a 
region  of  size  2(2Rq+D)  <<2 irc/w  is  that  corresponding  to  the  dipole  induced 
in  that  region  can  be  justified  rigorously. The  scattering  due  to  the 
planar  interface  is  treated  here  exactly  by  using  Fresnel  formulae  for  . 


r  i2i 

(3)  It  is  well  known1  1  that  the  magnetic  dipole  of  a  perfectly 
conducting  sphere  makes  an  important  contribution  to  the  intensity  of 
the  light  scattered  by  such  an  object.  The  electrostatic  calculation 
carried  out  by  us  eliminates  this  effect.  Fortunately,the  magnetic  dipole 
effects  are  important  in  the  radiation  zone  only,  and  contribute  negligibly 
to  the  near  field  computed  here. 


III.  The  results  of  the  numerical  calculations 


III.  1.  Objectives. 

The  electromagnetic  properties  of  the  structure  considered  here 
depend  on  a  large  number  of  parameters.  To  diminish  the  amount  of 
interesting-but-not-essential  information  we  plot  only  the  "enhancement 
ratio" 

t=  (E  E*  +E  E*  +  E  /  E  •  £  *  ,  (32) 

'  H  H  T)  V  <P  o  o  ’  '  ' 

which  is  the  ratio  of  the  local  and  the  incident  intensity.  The  graphs  pre¬ 
sented  in  what  follows  have  the  following  objectives. 

(1)  We  want  to  study  the  frequency  dependence  of  1  in  order  to 
learn  the  resonance  frequencies,  and  the  width  and  magnitude  of  the  resonant 
enhancement.  By  varying  the  geometrical  parameters  (i.e.  A.  =  D/Rq)  and 
nature  of  the  active  material  we  want  to  exemplify  the  "tunability"  of  the 
structure. 

(2)  In  analyzing  the  explicit  formulae  for  the  expansion  coefficients 

Am  and  Bm  we  discuss  the  origin  of  the  observed  resonances  and  their 
n  n 

frequency,  and  also  elucidate  bhetr  foie  Ly\  ^ivtY\Cj,  yu^e  to  observ  ed 

resonant  enhancement. 

(3)  Finally,  a  group  of  graphs  illustrates  the  spatial  properties  of 
the  enhancement  factor  I. 

III. 2.  The  materials. 

The  calculations  are  valid  for  systems  chosen  so  that  the  sphere  is 
a  good  conductor  and  its  size  satisfies  the  requirements  discussed  in  Section  II. 
There  are  no  limitations  on  the  active  materials  except,  of  course,  for  the 
fact  that  we  desire  large  enhancements  in  a  spectroscopically  interesting 
frequency  range;  this  is  accomplished  if  the  electromagnetic  resonances  of 
the  structure  that  appear  in  the  desirable  frequency  range  are  poorly  damped. 


Silver  is  by  now  the  best  loved  enhancer  for  incident  visible  (blue 
and  near  UV)  light.  The  flat  surface  has  a  surface  plasmon  with  the 
dispersion  relation  =  (co/c )  Re  ^€j(co)/fej(co)+l^.  This  cannot  be  excited 
by  light  but  it  cam  be  excited  by  any  localized  source  placed  near  the  surface, 
such  as  an  excited  molecule^*^  or  a  polarized  sphere. ^*^  Such  sources 
excite  the  high  k-  region  of  the  dispersion  curve  (i.e.  the  unretarded 

»  r 

plasmon)  at  the  frequency  given  by  Re  f^(co)=-l,  For  silver1 

this  is  at  ui  =3.62  eV. 
sp 

0- Silicon  carbide  is  a  newcomer  which  is  an  excellent  enhancer  in 

r  1 7] 

infrared.  Its  dielectric  constant  is  given  by1 

f  (to)  =  <00 -  to2  -  i toy  )  (tt£p-  to2  -  i coy  )”*  (33) 

The  parameters  appearing  in  this  equation  are  =6.7,  co^,  =  794  cm”*, 
(i.e.  the  restrahl  frequency  corresponding  to  maximum  reflectivity), 
co^  =  970  cm”*  and  y  -  8.5  cm”*.  The  surface  resonance  is  a  surface 
polariton.  Its  high  k..  (i.e.  unretarded)  frequency  is  given  by  e  (co)=-l 
and  equals  cogp  =  949  cm”  . 

InSb  is  a  semi-conductor  whose  dielectric  constant  is  given  by^*®^ 

Sto2  «_  e 

C{w)  =  +  j,  ■■  -  -  <34> 

coQ  -co  -i coy  co  +  xco  coc 

Hie  parameters  are:  c  =  15.7,  S  =  2,2,  coQ  =  179  cm"*,  y  =  2.7  cm  *, 

co  =  1023  cm”*,  oo  =30  cm"*, 
p  c 


The  second  term  in  Eq.  (34)  represents  contributions  to  the  dielectric 
constant  caused  by  optical  phonons.  The  third  is  due  to  conduction  electrons. 


A  useful  feature  of  n-type  InSbis  that  we  can  modify  the  bulk  plasmon  fre¬ 
quency  by  doping,  since  tOp  depends  on  the  free  electron  density  N  through 
=  (4irNe^/m*eQ0)^^. 

The  material  has  three  unretarded  surface  modes,  since  there 

are  three  frequencies  at  which  Re  e{o>)  equals  -1.  The  low  frequency 

mode  is  mainly  a  surface  polariton  and  the  high  frequency  one  is  mainly 

a  surface  plasmon.  Since  the  width  of  the  resonance  equals  ImC(co) 

taken  at  the  resonance  frequency,  only  the  plasmon  like  mode  is  expected  to 

play  an  important  role  if  high  local  fields  are  desired.  The  other  two  can, 

however,  be  very  effective  in  dissipating  the  energy  of  vibrationally  excited 

molecules  which  resonate  with  them. 

Doping  affects  the  plasmon  mode  by  shifting  its  frequency  to  higher 

1/2  - 1 

values  until  the  maximum  frequency  00  [ € / (c +1 )]  =  992  cm”  is  reached. 

p  L  03  00  J 

Through  plasmon-polariton  "interaction"  doping  also  shifts  the  frequency  of 
the  polariton  resonance.  The  minimum  value  is  U)Q. 

We  also  consider  an  "organic  material"  of  the  kind  studied  by  Pockrand 
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et.al.1  J  in  their  ATR  work.  A  representative  dielectric  constant  is 

2  f 

*<«>■  J  “  T-.-7-  •  (35> 

CdT-CO  ”1  CO/  0 

with  t  =  2.25,  f  =  0.2,  =2.24  eV  and  <5=  22.34  eV  The  unretarded 

surface  mode  is  a  surface  exciton. 


i 


III. 3.  The  Frequency  dependence  of  the  enhancement  factor. 


in. 3. a.  Silver  substrate. 

In  Fig.  2  we  plot  I,  as  a  function  of  frequency,  for  two  different 
values  of  the  geometrical  factor  A=D0/R.  The  point  at  which  T”  is 
calculated  is  located  on  the  flat  interface  directly  below  the  sphere 
(at  xo  =  zo  =  0  in  Fig.  1).  This  point  was  chosen  because  the  enhancement 
factor  there  is  very  high.  The  frequency  dependence  of  T  (i.e.  the  shape 
of  the  curve)  is  weakly  dependent  on  the  observation  point  as  can  be  seen 
from  the  Eqs.  (11)  and  (16)  —  (19).  The  observation  point  enters  in  the 
spherical  harmonics  and  the  exponentials  exp  [±  (n+l/2)/u].  The  resonant 
behaviour  (i.e.,  the  shape  of  the  frequency  dependence  curve)  enters  in 
A^1  and  B^1 .  A  change  in  observation  point  can  change  radically  the 
magnitude  of  the  enhancement  but  not  the  resonant  structure  of  the  curve 
(as  long  as  the  point  is  close  enough  to  the  sphere).  The  maximum  enhance 
ment  is  T  =  8300  (for  \=  0.05  given,  for  example,  by  RQ  =  100  A  and 
D  =  5  A)  and  it  occurs  at  3.27  eV.  The  resonance  is  very  broad.  If  the 
sphere  is  moved  further  from  the  flat  surface  the  maximum  enhancement 
is  (6.315)*  10  at  frequency  3.43  eV.  Note  that  in  the  case  of  Raman 
scattering  the  enhancement  of  the  emission  is  comparable  to  that  of  the 
local  field^^  and  therefore  if  I  ~10^  the  Raman  enhancement  is 


111.3.b.  £ -Silicon  carbide  substrate. 

•"W 

The  dependence  of  I  on  u>,  when  the  active  medium  ia  £-SiC,is 
displayed  in  Fig.  3a.  For  X=0.1  (i.e.  Rq  =  100  Ji ,  D  =  10  A  )  I  has  a 
peak  value  of  ~104,  located  it  u  :  928  cm’1.  For  0.05  (i.e  Ro*100A 
and  D  *  sk)  there  are  two  peaks,  at  oo  =  916  cm’1  and  935  cm’1.  This 
happens  because  the  damping  of  the  gap  modes  is  too  small  to  make  them 
merge  in  a  broad  peak.  This  is  further  discussed  in  Section  III. 4. 

The  s -polarized  ’ight  is  not  as  efficient  in  creating  large  local  fields. 

This  is  seen  in  Fig.  3(b),  where  I  is  plotted  for  conditions  identical  to 

Fig.  3(a),  except  for  the  polarization.  The  s-polarized  light  yields  an 

enhancement  of  '“600,  which  is  about  two  orders  of  magnitude  lower  than 

that  obtained  with  p-polarized  light.  This  is  the  result  of  a  symmetry 

effect  which  appears  because  the  electric  field  parallel  to  the  flat  surface 

does  not  excite  modes  in  which  the  induced  charge  oscillates  perpendicular 

to  the  surface.  Mathematically  this  appears  in  the  equations  through  the  fact 

that  A1”  =  Bm  =  0  for  m^  t  1  so  that  the  mode  m=0  is  not  excited.  The 
q  n 

p-polarized  light  excites  all  the  modes,  and  it  excites  the  m=0  mode  more 
strongly  as  9 —  90°.  Thus,  p-polarized  light  at  grazing  incidence  is  most 
effective. 

We  have  found  previously^ a  propensity  rule  which  states  that  a 
maximum  local  field  is  created  between  two  boundaries  if  the  electric  field 
of  the  excitation  source  oscillates  along  the  line  uniting  the  two  boundaries. 
This  happens  as  if  the  electrons  try  to  follow  the  field  and  are  pushed  against 
the  boundaries  causing  a  charge  concentration  which  increases  the  field  in  the 
vacuum,  between  the  boundaries. 


IH.3.C.  InSb  substrate 


The  frequency  dependence  of  I,  for  the  case  when  the  active  medium 
is  3hSb,  is  presented  in  Figs.  4(a)  and  (b).  The  incident  laser  character¬ 
istics  and  the  geometry  of  the  system  are  the  same  as  in  Fig.  2.  The 
figures  differ  through  the  bulk  plasmon  frequency  which  is  u)  =  195  cm"\ 
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for  '4(a^and  u)p  =  1023  cm”  ,for  4(b).  Such  differences  can  be  achieved 
through  doping.  At  low  conduction  electron  concentration  one  can  see  two 
sets  of  gap  modes,  one  at  140  cm  *  and  the  other  at  200  cm  These 
bands  are  composed  of  several  gap  modes  which  have  a  sizeable  width  and 
overlap.  The  shoulder  at  170  cm  *  is  an  indication  of  the  existence  of  such 
modes.  In  Fig.  4b  we  show  the  local  enhancement  factor  caused  by  the 
surface  plasmon  in  heavily  doped  LaSb.  The  surface  polariton  does  not 
give  a  large  local  enhancement,  probably  because  of  large  values  of  ImXj(co) 
at  that  frequency.  Note  that  doping  shifts  upwards  the  frequency  of  the  gap 
mode  and  increases  substantially  the  enhancement  at  that  frequency. 


IH.3.d.  Organic  substrate 

We  have  carried  out  a  calculation  with  the  active  medium  whose 

fig] 

dielectric  constant1  J  is  typical  of  organic  crystals  such  as  anthracene. 

The  surface  mode  excited  by  the  sphere  to  form  the  gap  modes  is  a  surface 
exciton.  The  enhancement  is  modest  (Fig.  5  )  as  compared  to  the  other 
active  media  considered  here,  but  it  is  sufficiently  large  to  carry  out  surface 
enhanced  spectroscopy. 


III. 4.  The  frequency  of  the  gap  modes  and  the  structure 
of  the  resonance  spectrum. 


In  this  section  we  wish  to  examine  more  closely  the  connection 

between  the  gap  mode  frequencies  (i.e.,  the  solutions  of  Eq.  (27))  and 

the  overall  frequency  response  of  the  system.  As  indicated  earlier,  the 

gap  mode  frequencies  are  determined  by  the  vanishing  of  the  denominators 

of  the  expansion  coefficients  A  and  B  .  These  denominators  are  of 

n  n 

the  form  |exp[(2n+l)^]  +  X1(o>)|,  where  is  given  by  Eq.(  9),  ^(oo)  = 

=  (1  -  CjM  (1  +  (to))  1  and  n  =  0,  1,  2. . .  If  we  plot  -Re\,  (co)  and 

-ImXj(co)  as  functions  of  the  (gap  mode)  resonance  frequencies  are 
located  at  the  intercept  of  the  horizontal  lines  corresponding  to  exp((2n  +  l)^) 
with  -ReXj(ui«  The  enhancement  for  each  resonance  is  inversely  pro¬ 
portional  to  the  value  of  ImXj(oj)  at  its  frequency.  This  is  illustrated  in 
Fig.  6  for  a  SiC  substrate.  The  dotted  lines  are  the  magnitudes  of 
exp  ((2n  +  l)^)  for  n  =  0,  1,  2  and  3.  For  n  =  0  there  are  two  resonances, 

one  at  cu  ^  833  cm-1  and  the  other  at  u:  ^  =“  949  cm-1  .  The  upper  mode 

o  o 

(2) 

tuo  occurs  at  practically  the  same  frequency  as  the  surface  polariton  of 
the  flat  substrate,  i.e.,  it  is  obtainable  to  good  accuracy  as  the  root  of  the 
equation  ReXj(co)  =  0  or  C^oj)  =  -!•  However,  unlike  the  latter,  it  is 
strongly  damped  (due  to  the  large  value  of  ImXj(a/^  and  is  therefore  not  a 
well  defined  excitation  of  the  sphere-plane  system.  The  lower  mode 
has  a  very  small  damping  associated  with  it  and  is  therefore  a  well  defined 
excitation  of  the  system.  The  role  of  the  conducting  particle  atop  the 
substrate  is  thus  seen  to  be  the  following:  it  strongly  damps  the  original 
(surface  polariton)  mode  at  and  creates  a  new  (gap-mode)  resonance 

at  the  appreciably  lower  frequency  These  features  persist  for  the 

higher  gap  modes.  For  n  =  1,  2  the  gap  modes  occur  at  929  cm"*  and 


940  cm  *,  respectively.  For  n  ^  3,  there  is  no  intercept  between 
-ReXj(a;)  and  exp  ((2n+l)^)  and  thus  there  are  no  well  defined  gap  modes. 
Another  way  of  saying  this  is  that  for  n  a  3  the  gap  modes  have  been  pushed 
up  towards  the  surface  polariton  frequency  of  949  cm  *  where  they  are 
strongly  damped  because  of  the  large  value  of  ImXj(ou). 

The  foregoing  analysis  would  seem  to  indicate  that  of  all  the  gap 
modes  it  is  only  the  lowest  few  (perhaps  the  first  two  or  three)  that  are 
important  in  determining  the  response  of  the  system.  In  fact,  in  the  case 
under  study  here,  the  lowest  (n=0)  mode  alone  seems  to  be  all  important 
since  its  damping  is  more  than  an  order  of  magnitude  smaller  than  for  the 
higher  modes.  One  might  therefore  be  tempted  to  conclude  that  the  overall 
resonant  response  of  the  system  would  be  at  the  frequency  of  the  lowest 
gap  mode  (=388  cm  *)  and  have  a  magnitude  inversely  proportional  to  the 
square  of  ImX^cc)  there.  The  actual  situation  turns  out  to  be  more  compli¬ 
cated  than  this.  The  reason  is  that  the  relative  contribution  of  the  different 
modes  to  the  potential  (and  hence  to  the  field  intensity)  is  given  by  an 
expression  of  the  form  (Eq.  (11)) 

nmA?exp  [(n  +  1/2W 

In  the  preceeding  analysis  we  looked  only  at  the  coefficients  A°. 
However,  in  order  to  obtain  the  relative  importance  of  different  modes  we 
must  compare  the  quantities  A^1  e^n  +  y1^  for  different  n.  This  is 

done  in  Table  1,  where  we  restrict  our  attention  to  terms  with  m  =  0.  In 
the  second  column  we  list  the  absolute  value  (modulus)  of  A°  .  As  expected, 
A°  decreases  rapidly  with  increasing  n  (the  frequency  at  which  these 
calculations  have  been  made  is  that  of  the  lowest  gap  mode  =  888  cm-*). 


In  the  3rd  column  we  show  the  quantity  e(n  +  l/2)Mo  and  in  the  4th  the 
product  e(n+*/2)Mo  |Ao  |  ^we  (jisregar^  here  the  spherical  harmonic 
Y°  which  is  generally  o^  order  unity).  Note  that  the  4th  column  decreases 
much  less  rapidly  with  increasing  n  than  the  2nd  column.  This  is 
because  the  decrease  of  |  is  offset  to  some  extent  by  the  increase  of 
the  exponential  factor  e^n  +  ^^^o,  The  4th  column  shows  clearly  that 
it  is  not  just  the  lowest  mode  which  is  important  but  a  large  number  of  low 
lying  modes;  this  is  evident  from  the  fact  that  the  cumulative  effect  of 
modes  1  through  7  is  as  important  as  the  lowest  (n=0)  mode.  For  n  larger 
than  10  or  so  |A°  |  plunges  sharply  and  so  does  the  product  in  the  last 
column,  thus  making  the  contribution  of  these  modes  rapidly  less  important. 

We  can  conclude  from  our  analysis  of  the  previous  paragraph  that 
the  overall  frequency  response  of  the  system  is  determined  by  the  contribution 
of  a  large  number  of  low-lying  gap  modes,  and  not  just  the  lowest  mode  as 
our  preliminary  analysis  seemed  to  indicate.  The  competition  between  the 
modes  results  from  the  finite  damping  of  the  substrate.  If  the  substrate  had 
no  damping  (Ime^uj)  =  0)  there  would  be  an  infinite  number  of  gap  modes  at 
each  of  which  the  response  (measured  by  the  local  field  intensity)  would  be 
infinite.  With  finite  damping  {Im€j(u;)  ^  0)  the  modes  are  broadened  and 
overlap  with  each  other.  As  a  result  of  this  smearing  out,  the  system 
response  peaks  at  a  frequency  intermediate  between  that  of  the  lowest  gap 
mode  (888  cm  *)  and  the  upper  surface  polariton  cutoff  (949  cm  *)  —  in  the 
present  case  this  is  at  916  cm  *  with  a  secondary  peak  at  935  cm”*. 


III. 5 .  The  spatial  extent  of  the  enhanced  field. 


We  present  now  several  graphs  which  display  the  spatial  distribution  of 

<N/ 

the  localized  enhancement  I  ,  when  the  excitation  frequency  corresponds  to 
the  maximum  local  field  enhancement.  We  take  \  =  0.05  (i.e.  Rq  =  100a 
and  D  =  5  A  ),  0Q  =  60°,  and  p-polarized  light.  We  compute  I  in  the  plane 
of  incidence  of  light,  along  specific  curves.  In  Figs.  8  we  show  the  variation 
of  I  along  a  circle  surrounding  the  sphere.  The  distance  to  the  sphere 
surface  is  1  A  •  The  variable  is  the  polar  angle  9  of  the  observation  point. 
The  largest  enhancement  is  obtained  for  9  =  -  180°,  which  corresponds  to 
points  located  between  the  sphere  and  the  planar  surface.  The  field  has  a 
minimum  near  9=  -90°,  which  is  located  in  the  "shadow"  of  the  sphere. 
Larger  fields  are  obtained  on  the  "sunny"  side  of  the  structure  (i.e. 

0  <  9  1  90). 

Note  that  the  spatial  distribution  of  the  local  intensity  is  very  similar 
for  Ag,  SiC  (Fig.  8(a)),  IhSb  and  the  organic  material  (Fig.  8(b)).  All  these 
graphs  were  computed  at  the  resonance  frequency  of  each  material.  The 
curves  would  probably  be  very  different  if  they  were  all  computed  at  the  same 
frequency. 

In  Fig.  9  we  show  the  fields  along  the  symmetry  axis  of  the  system, 
in  the  space  between  the  sphere  and  the  plane.  The  enhancement  is  very 
high  and  constant  along  the  line.  In  Fig,  10  we  plot  the  enhancement  factor 
along  the  line  AB  shown  in  Fig.  1.  The  local  intensity  decays  as  one  moves 
from  A  towards  B.  At  points  approximately  100  a.  away  from  A,  I  reaches 
the  small  value  it  would  have  for  a  flat  surface  in  the  absence  of  the  sphere. 
The  curves  are  very  similar  for  different  active  materials  as  long  as  the 
calculation  is  carried  out  at  the  resonance  frequency  of  each  material. 


IV .  Discus  sion 


The  calculations  presented  here  indicate  that  the  gap  modes  formed 
between  a  perfectly  conducting  sphere  and  an  optically  active  medium  are 
potentially  useful  for  surface  enhanced  spectroscopy.  The  sphere  permits 
the  optical  excitation  of  the  flat  surface  mode  (surface  plasmon,  or  polariton, 
or  exciton)  but  it  distorts  it  so  strongly  that  it  seems  useful  to  regard  these 
resonances  as  new  modes.  The  name  of  gap  modes  seems  appropriate.  The 
local  intensity  is  concentrated  between  the  sphere  and  the  plane,  with  large 
contributions  at  the  illuminated  side  of  the  sphere  and  smaller  contributions 
on  the  shady  side.  This  is  very  different  from  a  mode  of  the  flat  surface 
which  is  delocalized  along  the  plane. 

The  enhancement  factor  produced  by  the  gap  modes  is  very  large, 
but  unfortunately  its  spatial  extent  is  small.  This  is  not  by  accident,  since 
one  can  think  of  the  enhancement  process  as  resulting  from  trapping  a  photon 
in  the  gap  region.  Thus  the  delocalized  photon  energy  is  packed  in  a  small 
volume  and,  since  the  field  energy  per  unit  volume  is  proportional  to  E  •  E  ,  the 
conservation  of  the  total  energy  re.quires  that  E  •  E  be  large  in  the  small  region. 
The  ability  of  generating  large  enhancements  is  insured  by  large  oscillator 
strength  (so  that  the  number  of  absorbed  photons  is  large)  and  the  existence 
of  electromagnetic  and  geometric  factors  which  pack  the  photon  energy  in  a 
small  space.  Our  structure  seems  to  satisfy  these  requirements.  However, 
the  number  of  molecules  that  can  be  exposed  to  high  fields  goes  down.  One 
can  try  to  compensate  for  this  by  depositing  on  the  surface  a  large  number  of 
spheres,  with  the  largest  radius  for  which  ZRo«2rr  c/ui.  However,  one 
should  keep  in  mind  that  if  the  sphere  concentration  becomes  too  large  the  gap 
modes  may  interact  with  each  other,  become  delocalized  and  thus  lower  the 
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enhancement  factor.  Such  effects  are  very  difficult  to  compute  accurately  or 
estimate  reliably;  it  is  probably  better  to  study  them  experimentally. 


One  can  envision  several  types  of  applications,  (a)  Pre-adsorb  the 
molecules  of  interest  on  an  active  single  crystal  substrate  (e.g.  Ag,  /3-SiC, 
InSb)  and  deposit  on  the  flat  surface  chemically  inert  perfect  metal  spheres, 
imbedded  in  an  optically  and  chemically  inert  matrix.  One  can  then 
attempt  to  use  the  structure  thus  obtained  to  study  the  spectroscopic  properties 
of  the  preadsorbed  molecules,  (b)  Assume  now  that  one  is  interested  in  the 
spectroscopic  properties  of  molecules  adsorbed  on  small  metal  particles. 

To  carry  out  enhanced  IR  spectroscopy  one  can  place  these  particles  on  a  InSb 
surface  alter  adsorbing  on  them  the  molecules  of  interest.  One  can  hope,  by 
doping  InSb,  to  be  able  to  tune  the  structure  to  have  a  resonance  at  the 
molecular  frequency  of  interest.  The  molecules  will  absorb  very  strongly  at 
the  molecular  frequency  and  thus  depress  the  resonance.  One  can  hope  to 
detect  this  effect  in  several  ways:  (a)  difference  spectra,  with  and  without 
the  absorbing  molecules:  (0)  polarization  modulation,  since  the  absorption 
from  a  p-polarized  beam  is  much  stronger  than  that  from  a  s -polarized 
beam;  (y)  by  monitoring  the  local  temperature  with  a  bolometer.  We  are 
now  in  process  of  carrying  out  computations  to  evaluate  the  chances  that  such 
measurements  might  succeed,  (c)  The  structure  can  also  be  used  to  attempt 
to  carry  out  infra-red  (on  InSb  and  /3-SiC)  or  visible  (Ag)  multiphoton  photo¬ 
chemistry  or  non-linear  spectroscopy,  (d)  Finally,  one  can  think  of  making 
use  of  the  effects  described  here  to  construct  IR  detectors,  especially  if 
one  is  interested  in  having  a  detector  for  a  small,  given  frequency  range. 

One  can  use  an  InSb  film  and  metal  particles  trapped  on  its  surface  by  a 
matrix  whose  molecules  absorb  in  the  desired  range.  By  doping  the  InSb 


film  one  can  tune  the  gap  modes  to  the  desired  frequency.  One  can  then 
detect  ER  radiation  of  the  specified  frequency  by  monitoring  the  film 
temperature. 

We  emphasize  that  these  suggestions  are  quite  speculative  and  that 
we  mention  them  only  because  they  deserve  further  contemplation  and/or 
experimental  testing. 
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Figure  captions 


Fig.  1.  The  geometry  of  the  system:  A  planar  interface  separates  the 

active  medium,  characterized  by  a  complex  frequency- 

dependent  dielectric  constant  Cj(uj),  from  an  "inert"  medium 

of  dielectric  constant  e  •  Imbedded  in  the  inert  medium  at 

o 

distance  D  above  the  interface  is  a  perfectly  conducting 
spherical  particle  of  radius  Rfl  .  The  incoming  laser  beam 
has  a  wavevector  k  and  electric  field  Eq.  The  figure  shows 
a  p-polarized  beam. 


Fig.  2.  Resonance  spectrum  of  a  Ag  substrate  with  a  perfectly  conducting 
sphere  on  it.  The  intensity  enhancement  ratio  T  (defined  in  text) 
has  been  plotted  as  a  function  of  the  incident  laser  frequency  in 
eV.  The  laser  beam  is  incident  at  angle  Bq  =  45°  and  is 
p-polarized.  The  full  line  is  computed  for  \  =  0.05  (e.g., 

Ro  =  100  A  ,  D  =  5  A  );  the  dott  ed  line  is  for  \  =  0.1  (e.g., 

Ro  =  100 A,  D  =  10 A).  The  value  Of  FA  indicated  for  each 
curve  is  the  factor  by  which  the  vertical  axis  must  be  multiplied 
to  obtain  the  value  of  I* for  that  curve.  This  convention  is  followed 
on  all  subsequent  figures. 


Fig.  3a.  Same  as  Fig.  2,  but  the  substrate  is 

Fig.  3b.  Same  as  Fig.  2,  but  the  substrate  is 
light  is  s -polarized. 


0-  SiC 

0-SiC  and  the  incident 


Fig.  4a.  Same  as  Fig.  2,  but  the  substrate  is  n-type  InSb;  the  substrate 


is  lightly  doped  and  has  a  free  electron  concentration  corresponding 
to  a  plasma  frequency  of  cOp  =  195  cm  *. 

Fig.  4b.  Same  as  Fig.  4a,  but  with  a  heavily  doped  substrate  corresponding 

to  co  =  1023  cm” * 

P 

Fig.  5.  Same  as  Fig.  2,  but  with  an  organic  substrate. 

Fig.  6.  The  plot  of  ImX  (<o)  and  Re  X  (co)  for  a  SiC  substrate  for  X=  0.05. 

The  horizontal  lines  represent  exp  [(2n+l)/iQ]  for  n=0,  1,2,  3. 

See  Section  III. 4  for  a  further  discussion. 

Fig.  7.  Spatial  variation  of  enhancement  ratio  I  at  resonance. 

Log^T  is  plotted  as  a  function  of  the  angular  coordinate  9  for 

fixed  f  =  1A.  The  full  line  is  for  a  Ag  substrate  and  co  =  3.27  eV; 

dcLS^-iot  _1 

the  dufatesaefr  line  is  for  a  SiC  substrate,  co  =  916  cm  .  The 

system  parameters  are  X  =  0.05  (e.g.,  R0  =  100  A,  D  =  5  A). 

The  external  laser  is  incident  at  0Q  =  60°  and  is  p-polarized. 

Fig.  8.  Same  as  Fig.  7  but  with  the  following  substrates:  dotted  line 
for  a  n-type  InSb  (cOp  =  1023  cm"*),  and  co  =  931  cm  *;  the 
hashed  line  is  for  an  organic  material  substrate  and  co  =  2.27  eV. 

Fig.  9.  k°8io  1  plotted  as  a  function  of  the  coordinate  ZOfor  XO  =  0 
(along  z-axis).  The  system  parameters  and  substrates  used 
are  identical  to  those  in  Figs.  7  and  8. 

0 

Fig.  10.  Log  10  1  plotted  as  a  function  of  the  co-ordinate  XO  for  ZO  =  1A  . 
The  system  parameters  and  substrates  used  are  identical  to 
those  in  Figs.  7  and  8. 


Table  1. 

An  analysis  of  the  contributions  of  the  gap  modes 
to  the  frequency  dependence  of  the  local  intensity. 
See  Section  III. 4  for  discussion. 
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